How do you know if you ran through a wall? 
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Stable topological defects of light (pseudo) scalar fields can contribute to the Universe's dark 
energy and dark matter. Currently the combination of gravitational and cosmological constraints 
provides the best limits on such a possibility We take an example of domain walls generated by 
an axion-like field with a coupling to the spins of standard-model particles, and show that if the 
galactic environment contains a network of such walls, terrestrial experiments aimed at detection of 
wall-crossing events are realistic. In particular, a geographically separated but time-synchronized 
network of sensitive atomic magnetometers can detect a wall crossing and probe a range of model 
parameters currently unconstrained by astrophysical observations and gravitational experiments. 

PACS numbers: 14.80.Va, 98.80.Cq 



Introduction. Despite a remarkable success of the 
Standard Model in describing all phenomena in parti- 
cle physics, the cosmology presents a formidable puzzle, 
with dark energy and dark matter - two substances of un- 
known origin - comprising about 75% and 20% of the Uni- 
verse's energy budget. Last decades have seen a dramatic 
expansion of all experimental programs aimed at clari- 
fying the nature of dark matter (DM) and dark energy 
(DE). While many widely ranging theories of dark matter 
exist, most of the experimental efforts go into searches of 
dark matter of some particle physics variety, producing 
upper limits on the DM-atom interaction strength. Tests 
of DE models occur on cosmological scale, showing so far 
its consistency with the cosmological constant. 

The purpose of this Letter is to show that a new class 
of objects, stable topological defects (such as monopoles, 
cosmological strings or domain walls), that will con- 
tribute both to the DM and DE, can be searched for and 
studied with the global network of synchronized atomic 
magnetometers. To be more specific, we consider an ex- 
ample of a domain wall network built from the axion-like 
fields. Our focus on axion-like fields and the pseudoscalar 
interaction of these fields with matter, is motivated by 
the theoretical considerations of "technical naturalness" , 
that allow preserving the lightness of the pseudoscalar 
fields despite a significant strength of interaction with 
matter. Observable effects of light pseudoscalar particles 
can vary considerably, depending on their mass m a . We 
refer the reader to a sample of literature on the subject, 
covering a wide range of m a from 10 -33 to 10 5 eV pQ. 

Scalar-field potentials with some degree of discrete 
symmetries admit domain-wall-type solutions interpolat- 
ing between domains of different energy-degenerate vacua 
[2]. In these models, initial random distribution of the 



scalar field in the early Universe leads to the formation of 
domain- wall networks as the Universe expands and cools. 
For Q CD-type axions, if stable, such domain walls could 
lead to disastrous consequences in cosmology by storing 
too much energy [2]. For an arbitrary scalar field, where 
parameters of the potential are chosen by hand, the "dis- 
aster" can be turned into an advantage. Indeed, over 
the years there were several suggestions how a network 
of domain walls could be a viable candidate for DM or 

de ng. 

Herein, we revisit a subset of these ideas from a prag- 
matic point of view. We would like to address the follow- 
ing questions: (1) if a network of domain walls formed 
from axion-like fields exists in our galaxy, what are the 
chances for an encounter between the Solar system and a 
pseudoscalar domain wall? and (2) how could the event 
of a domain-wall crossing the Earth be experimentally 
determined? Given gravitational constraints on the av- 
erage energy density of such walls and constraints on the 
coupling of axion-like fields to matter [5 , 6 , it is not ob- 
vious that the allowed parameter range would enable de- 
tection. Yet we show in this Letter that there is a realistic 
chance for the detection of the domain walls, even when 
the gravitational and astrophysical constraints are taken 
into account. This goal can be achieved with correlated 
measurements from a network of optical magnetometers 
with sensitivities exceeding 1 pT/vTtz, placed in geo- 
graphically distinct locations and synchronized using the 
global positioning system (GPS). 

Physics of light pseudoscalar domain walls. We start 
by considering the Lagrangian of a complex scalar field 
(/>, invariant under Z^-symmetry, <\> — >• exp(z27rfc/7V)0, 
where k is an integer. We choose the potential in such a 
way that it has N distinct minima 
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where So has dimension of energy and A is dimensionless. 
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Choosing <p = 2 -1 / 2 S exp(ia/ So) to parameterize the 
scalar field, we find that the potential V{<j>) is minimized 
for the following values of S and a, 



S = S ; a = S x jo; 
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Freezing the Higgs mode to its minimum, S — So, pro- 
duces the effective Lagrangian for the a field, 
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with Vo = 4ASq. The spatial field configuration a(r) 
interpolating between two adjacent minima represents a 
domain-wall solution. A network of intersecting domain 
walls is possible for TV > 3. The solution for a domain 
wall along the xy-plane that interpolates between a = 
and 2ttSo/N neighboring vacua with the center of the 
wall at z = takes the following form, 
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The characteristic thickness of the wall d is determined 
by the mass m a of a small excitation of a around any 
minimum, d ~ 2/m a . The mass m a can be expressed 
in terms of the original parameters of the potential, 
m a = A^S " 1 (l/o/2) 1/2 = (2X) 1 / 2 NS . Owing to the fact 
that V(4>) can have many different realizations other than 
0, we shall use solution Q as an example, rather than 
a generic domain- wall profile for TV > 3. The important 
parameters are the gradient of the field inside the wall, 
m a So/TV, and m a , which determines the wall thickness. 

Gravitational and astrophysical constraints. From the 
macroscopic point of view at distance scales much larger 
than d, the wall can be characterized by its mass per 
area, referred to as tension, 
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The network of domain walls will have an additional 
distance-scale parameter L, an average distance between 
walls, or a characteristic size of a domain. This param- 
eter is impossible to calculate without making further 
assumptions about the mechanisms of wall formation 
and evolution. We treat it as a free variable and con- 
strain the maximum energy density of the domain walls, 
Pdw ~ cr/L in the neighborhood of the Solar System by 
the dark-matter energy density, pdm — 0.4 GeV/cm 3 , 



Pdw < Pbm 
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This constraint implies some flexible evolution of the 
domain- wall network and the possibility for them to build 
up their mass inside galaxies. We consider such the con- 
straint as the most conservative, i.e. giving the most 



relaxed bound on Pdw- If the network of domain walls 
is "stiff" and its density inside galaxies is not enhanced 
relative to an average cosmological value, then a stronger 
constraint can be derived by requiring that domain walls 
provide a (sub) dominant contribution to the dark-energy 
density, pow < Pde, where pde — 0.4 x 10 -5 GeV/cm 3 
[9]. In that case the constraint on So /TV is strength- 
ened by ~ 300. Our choice of the normalization for L 
and m a in (|6| is suggested by the requirement of hav- 
ing wall crossings within ~10 yr with relative velocity of 
v = 10 _3 c typical for galactic objects, and having the 
signal duration in excess of 1 ms. This choice can be 
self-consistent within the cosmological scenario for the 
formation of the domain- wall network from randomly dis- 
tributed initial a- ini assuming that the network is "frus- 
trated", and exhibits Pdw °° R 1 scaling, where R is the 
cosmological scale factor. As a word of caution, we add 
that the numerical simulations of domain walls in some 
scalar field theories have shown much faster redshifting of 
Pdw, and never achieved the frustrated state [7]. In light 
of this, some unorthodox cosmological/astrophysical sce- 
narios for the formation of domain walls may be required. 

We consider two types of pseudoscalar coupling of the 
field a with the axial- vector current of a standard-model 
fermion, J M = ^7^75^, 
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where the arrows show the reduction of these Lagrangians 
at the minima of V(a), and f- are free parameters of 
the model with dimension of energy. The normalization 
is chosen in a way to make connection with axion litera- 
ture. The derivative nature of these interactions softens 
problems with "radiative destabilization" of m a . It is 
also important that the effective energy parameters nor- 
malizing all higher dimensional interactions in ([7| and 
([8| are assumed to be above the weak scale. Both Cu n 
and £ qua d lead to the interaction of spins of atomic 
constituents and the gradient of the scalar field, 
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For light scalars of interest, the astrophysical bounds 
limit \f n ,p,e\ > 10 9 GeV [6], while bounds on quadratic 
d^a 2 interactions are significantly weaker, f[ > 10 TeV 
0. In what follows we will derive the signal from fa in 
d9j), and then generalize it to the j[ case. 

Spin signal during the wall crossing. The principles 
of sensitive atomic magnetometry are, for example, de- 
scribed in Ref. [10]. A typical device would use param- 
agnetic atomic species such as K, Cs, or Rb by them- 
selves or in combination with diamagnetic atoms whose 
magnetic moments are generated by nuclear spin (e.g., 
the spin-exchange-relaxation- free [SERF] 3 He-K magne- 
tometer described in Ref. [XT]). Specializing ([9| for the 
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case of two atomic species, 133 Cs in the F = 4 state and 
3 He in the F =1/2 state, we calculate the energy differ- 
ence AE between the F z = F and F z = —F states in the 
middle of the wall, 
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In these formulae we assumed that the nuclear spin is 
mostly due to unpaired neutron ( 3 He) or g 7 / 2 valence 
proton ( 133 Cs), and one can readily observe complemen- 
tary sensitivity to fi in two cases. We can express these 
results in terms of the equivalent "magnetic field" inside 
the wall using /iB e ffF/F = VaF / '(F '/eff) identification, 
where \i is the nuclear magnetic moment. The magnitude 
of B e R (direction is impossible to predict) is given by 
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and the larger equivalent field strength for 3 He originates 
from its smaller magnetic moment. The couplings and 
wall parameters in Eq. (11) are normalized to the max- 
imum allowed values from Eq. |6|. The duration of the 
signal is given by the ratio of wall thickness to the trans- 
verse component of the relative Earth-wall velocity, 
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Such crossing time can easily be in excess of the Cs mag- 
netometer response time t r , and we can combine the 
B™ a * and At into a signal factor S = B^(At) 1 ^ 2 to 
be directly compared to experimental sensitivity, 
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where in the inequality we used the gravitational con- 
straint from Eq. ^ . The maximally allowed value for the 
signal (rsj pT/Vnz), after taking into account the gravi- 
tational and astrophysical constraints, exceeds capabili- 
ties of modern magnetometers that can deliver FT/y/Hz 
sensitivity [10 . For the 3 He-K SERF magnetometer, the 
more appropriate figure of merit would be the tipping 
angle of the helium spin after the wall crossing, assum- 
ing that the typical crossing time is below the dynamical 
response time. Taking the spins to be oriented parallel 
to the wall, we calculate this angle to be 
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This could be far in excess of 10-nrad tipping angles 
that can be experimentally detected [12 . Thus, both 
types of magnetometers offer ample opportunities for a 
realistic detection of the wall-crossing events. So far we 
have used the galactic constraints (J6|, pow < Pdm- It 
is noteworthy that even if the energy density of walls 
in the galaxy does not exceed cosmological dark-energy 
density, i.e. Pdw < Pde, the expected signal can reach 
AO ~ 10" 5 rad and S - fT/VHz, which is still a realis- 
tic signal for detection with the best magnetometers. It 
is remarkable that a possible domain-wall component of 
DE can, in principle, be detected in the laboratory. 

Going over to /' couplings, we notice that the struc- 
ture of the signal is different: 5^ ax now changes sign, 
vanishing in the middle of the wall. Taking i^ ax at 
a = Sott/(2N) inside the wall, and skipping intermedi- 
ate states in a similar derivation, our sensitivity formulae 
( 11 ) and ( fl3| ) are modified according to the following sub- 
stitution, 
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where again /' is normalized on its minimum allowed 
value. One can observe a dramatic increase in the pos- 
sible signal due to a much weaker astrophysical con- 
straints on £ qua d. In Fig. 1, we plot the experimen- 
tal accessible parameter space in terms of characteris- 
tic time between wall crossing events, T = L/(10 _3 c), 
and strength of the coupling constants, / and /', fixing 
m a = 10 -9 eV,v±/c = 10 -3 for concreteness, and satu- 
rating either DM or DE density constraints. We assume 
that the magnetometer sensitivity is S = iT/a/Hz. The 
light (dark) shaded areas indicate the coupling range that 
can be realistically probed with the magnetometer net- 
work when DM(DE) constraints are saturated, by impos- 
ing all constraints and additionally requiring T < 10 yr. 
One can see that the large part of the parameter space 
is accessible, and for the case of C n ua,d even the DE con- 
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FIG. 1: Parameter space open for detection of the wall crossing, 
T/(lCT 3 c) in yr vs /(/' ) in GeV. Shaded triangle on the left cor- 
respond to £quad case, and on the right to Ain- 

Network of synchronized magnetometers. While a sin- 
gle magnetometer is sensitive enough to detect a domain- 
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wall crossing, due to the rarity of such events it would 
be exceedingly difficult to confidently distinguish a signal 
from false positives induced by occasional abrupt changes 
of magnetometer- operation conditions, e.g., magnetic- 
field spikes, laser-light-mode jumps, etc. A global net- 
work of synchronized optical magnetometers is an attrac- 
tive tool to search for galactic/cosmological domain walls, 
as it would allow for efficient vetoes of false domain- wall 
crossing events. 

Ideally, one would require n > 5 magnetometer sta- 
tions in such a network. The difference in timing ti of a 
putative signal is related to the transverse velocity and 
the unit normal vector to the wall, n, ti —tj = L>ij •nvj 1 , 
where L^- are the three-vectors of the relative positions 
of magnetometers i and j. Four stations are required to 
specify magnetometer-defined 3D system of coordinates, 
and three time intervals between four ti will enable to 
unambiguously determine the three- vector nvj 1 . This 
makes the predictions for the timing of the event at the 
fifth station, £5, which can be used as a tool for reject- 
ing accidental backgrounds. Consider a network of simi- 
lar magnetometers with fast response time separated by 
distances of 0(300 km) operating during a long period 
T ~ yr. Suppose that r is an average time between 
accidental spikes in the background above certain value 
B®ff that cannot be distinguished from the signal. Then 
the probability of having four events in four different sta- 
tions within time intervals corresponding to the typical 
wall travel time from station to station, t trav ~ l/v ~ s, 
is A234 ~ T£t rav r -4 , where we take T > r > ttrav 
To have this probability below one, one should achieve 
r > 100 s. If indeed four accidental background spikes 
lead to false signals in four stations within t tra v, the do- 
main wall interpretation will predict the event in the 
fifth station within a narrow window of the wall cross- 
ing At ~ ms, and the probability of this to happen due 
to accidental background is P12345 ~ (At/r)Pi234, or less 
than 10 -5 for r ~ 100 s. Increasing the number of sta- 
tions will enable to search for weaker signal and 
tolerate shorter r [15] . 

Recently we set up a prototype for the magnetome- 
ter network consisting of two magnetometers operated in 
magnetically shielded environments located in Krakow, 
Poland and Berkeley, USA (a separation distance of 
about 9000 km) . One of the magnetometers (Krakow) is 
based on nonlinear magneto-optical rotation [14], while 
the other magnetometer (Berkeley) is a SERF device 
[T5] . The magnetometers achieved comparable sensitivi- 
ties of 10 fT/Hz 1 / 2 , which can be further improved upon 
optimization. The expected parameters of the signal, 



At ~ 1 ms and the minimum time-separation between 
the events At tra v ~ 30 s, can be precisely determined us- 
ing a GPS time source (for more details see Ref. [16 ). We 
have recently performed proof-of-principle experiments 
[T6] demonstrating the ability to correlate the signals of 
two magnetometers. In particular, we demonstrated sig- 
nificant reduction of noise and rejection of false-positive 
events present in magnetometer signals. The measure- 




FIG. 2: Schematic of Global Network of Optical Magnetometers 
for Exotic physics (GNOME) planned to be used for detection of 
the domain-wall crossing. The wall-crossing events recorded with 
four magnetometers at ti allow determination of the normal veloc- 
ity of the wall v±. The remaining magnetometer(s) will be used to 
verified the measurements by predicting the time of the events in 
the locations (see text). 



ments proved the feasibility of correlated magnetic-field 
measurements opening avenues for further investigations 
involving more magnetometers. 

Summary. We have shown that a network of mod- 
ern magnetometers offers a realistic chance for detecting 
the event of an axion-type domain-wall crossing and can 
probe parts of the parameter space where such walls can 
contribute significantly to dark matter/dark energy. 
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